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Abstract 

We show that the combined data from solar, long-baseline and reactor neutrino experiments 
can exclude the generalized bicycle model of Lorentz noninvariant direction-dependent and/or 
direction-independent oscillations of massless neutrinos. This model has five parameters, which 
is more than is needed in standard oscillation phenomenology with neutrino masses. Solar 
data alone are sufficient to exclude the pure direction-dependent case. The combination of 
solar and long-baseline data rules out the pure direction-independent case. With the addition 
of KamLAND data, a mixture of direction-dependent and direction-independent terms in the 
effective Hamiltonian is also excluded. 
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1 Introduction 



Neutrino oscillations are now a well-established phenomenon [lj. Data from solar, atmospheric, 
reactor and accelerator experiments may be explained by the now standard scenario with three 
active, massive neutrinos, with the possible exception of the LSND experiment [2j. Recently it has 
been suggested that Lorentz-invariance and CPT violating interactions originating at the Planck 
scale can also lead to neutrino oscillations, with or without neutrino mass [31 HI El [7] . These 
interactions can be nonisotropic, which could lead to different oscillation parameters for neutrinos 
propagating in different directions. 

The effective hamiltonian that describes the evolution of massless neutrinos in vacuum in the 
presence of Lorentz-invariance violating interactions may be written as [5] 



where p^ = (E, —Ep) is the neutrino four-momentum, p the neutrino direction, and i,j are flavor 
indices. The coefficients ai have dimensions of energy and the cl are dimensionless. The Kronecker 
delta term on the right-hand side of Eq. (pQ) may be ignored since oscillations are insensitive to terms 
in h e ff proportional to the identity. For antineutrinos, ai — ► — cll- 

Direction dependence of the neutrino evolution enters via the space components of and cl- 
The coefficients may be specified in a celestial equatorial frame (T,X,Y,Z), which has Z axis 
along the Earth's rotation axis and X axis towards the vernal equinox. The two-parameter bicycle 
model [5] can be defined as follows: cl is isotropic, with only one nonzero element (cl)^ t = 2c, 
and (oi)e M = (ai)er = (0, aZ/y2) are the only nonzero a^. We have generalized the model by 
letting (ai)e A = (a^er = (0, an/v2), where h is the preferred direction for the interaction. This 
increases the number of parameters in the model to four, which is equal to the number required 
in the usual massive neutrino description of oscillations (two mass-squared differences and two 
mixing angles) pp. We also consider a five-parameter model which has a linear combination of 
direction-dependent and direction- independent a^. 

In this letter we examine the phenomenology of this direction dependence in the generalized 
bicycle model with massless neutrinos. We find that the pure direction-dependent bicycle model 
is ruled out by solar neutrino data alone, while a combination of solar and long-baseline neutrino 
data excludes the pure direction-independent case. A mixture of direction-dependent and direction- 
independent terms is excluded when KamLAND data are added. In Sec. 2 we present the model 
and the neutrino oscillation probabilities. In Sec. 3 we discuss the constraints from atmospheric 
and long-baseline neutrino experiments. In Sec. 4 we discuss the constraints from solar neutrino 
experiments, and in Sec. 5 we discuss the combined constraints, including KamLAND. In Sec. 6 we 
present our conclusions. 
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2 Neutrino oscillations in the generalized bicycle model 

Neutrino oscillations occur due to eigen energy differences in h e ff and the fact that the neutrino 
flavor eigenstates are not eigenstates of h e ff. For massless neutrinos p^ = (E, —Ep), where p is the 
direction of neutrino propagation. Then for the generalized bicycle model 
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where 



cos Q = p ■ h , (3) 
i.e., is the angle between the neutrino momentum and the preferred direction. From the diago- 



nalization of h e ff, there are two independent eigenenergy differences Ajk = Ej — Ej., 
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A 2 i = ^ ^EZ + Eftco&Q + E^ , A 32 = ^ ^E 2 + E 2 cos 2 8 - E^j , 



where rnk and Eq are defined in terms of the Lorentz-invariance violating parameters by 
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and the energy-dependent mixing angle is 
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The off-diagonal oscillation probabilities are [5] 



2 sin 2 9 cos 2 9 sin 2 (A 3 iL/2) 
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P(vp Vt) = sin 2 0sin 2 (A 2 iL/2) - sin 2 9 cos 2 9 sin 2 (A 3 iL/2) + cos 2 9 sin 2 (A 32 L/2) 

where A 3i = A 32 + A 21 . 

If Eq <C E 2 , i.e., a 2 <C (cE) 2 , for atmospheric and long-baseline neutrinos, then sin 2 9 <C 1, 
cos 2 9 ~ 1 and the only appreciable oscillation is 



where 



P(z^~^)~sin 2 (A 32 L/2) 



A 32 ~ J§ cos 2 6 . 



Thus the oscillation amplitude is maximal, the effective mass-squared difference is 

<5m 2 jy = itiq cos 2 , 



(9) 
(10) 
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and the energy dependence in this limit is the same as for conventional neutrino oscillations due to 
neutrino mass differences. Since the measured values for Sm^f agree for atmospheric neutrinos and 
the K2K [S] and MINOS [9j long-baseline experiments, the effective cos 2 G must also have similar 
values in all of these experiments. 

3 Atmospheric and long-baseline neutrinos 

3.1 Directional dependence 

With the Earth's rotation axis chosen as the Z direction and the position of the detector given by 
(0, (f>) in a standard spherical polar coordinate system (see Fig. [T]), the neutrino direction can be 
written as 

p = — cos r + sin (3{— sin a 9 + cos a (p) , (12) 

where r denotes the detector position, and the unit vectors f, 9 and <j> point in the upward, southerly 
and easterly directions, respectively. The angle (3 is the usual zenith angle (/? = for a downward 
event) and a denotes the compass direction of the neutrino velocity projected on the plane tangent 
to the Earth's surface (a = for a neutrino going in the eastward direction). We take the preferred 
direction to be 

n = sin £ cos % X + sin £ sin x Y + cos £ Z . (13) 
In our spherical polar coordinate system 

n = [sin£cos(0 + x) cos#l + cos£ sin#£,] f + [sin£ cos(0 + x) sin#£ — cos£cos#l] 9 

— sin £ sin(</> + %) <l> > (14) 

where the usual angle spherical polar 9 has been replaced by the latitude of the detector 9l = \ — 9. 
(positive for the northern hemisphere, negative for the southern hemisphere). The azimuthal angle 
<j) is chosen so that (ft = corresponds to the preferred direction so that the angle x may be 
dropped. The angular dependence in the oscillation formulas is then 

cos0 = cos <!; (sin /3 sin a cos #l — cos /3 sin 0^) 

— sin£ cos 0(sin /3 sin a sin 0l + cos (3 cos 9l) — sin£ sin (3 cos a sin^ . (15) 

In Eq. (115p . £ gives the orientation of the preferred axis with respect to the Earth's rotation axis, 
a (compass direction) and (3 (zenith angle) relate to the neutrino direction, and (j) depends on the 
time of the sidereal day (<p = when the detector is facing the preferred direction). 

To help understand the complicated angular dependences in Eq. (fT5|) . we consider three special 
cases: 

downward(/3 = 0) : cos = — (cos^ sin#£ + sin£ cos 9l cos <b) , (16) 
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A North 




Figure 1: Detector position in atmospheric and long-baseline experiments. The angle Ol is the 
latitude, while (f) (not shown) measures the time of the sidereal day. The unit vectors f,9,cp define 
the upward, southerly and easterly directions, respectively, for a neutrino event in the detector. 

upward(/3 = n) : cos O = cos£ sin6>£ + sin£ cos Ol cos 4> , (17) 
horizontal(/? = vr/2) : cos = cos£cos#£ sina — sin£(sin#c cos (j>sma + sin^cos a) . (18) 

Note that since only cos 2 appears in the oscillation formulas, the oscillation wavelengths for 
upward and downward events are the same. 

3.2 £ = 

If the preferred direction is aligned with the Earth's rotation axis, then £ = and 

cos 2 = (sin j3 sin a cos Ol — cos j3 sinOi) 2 . (19) 



Note that in this case does not depend on time of day (measured by (f)). For accelerator experi- 
ments with relatively short baselines compared to the Earth's radius (such as K2K and MINOS), 
the neutrino path can be considered to be in the plane that is tangent to the Earth's surface, so 
that Eq. (|18|) applies and cos 2 = sin 2 a cos 2 6l- Since the direction of the neutrino path in K2K is 
approximately given by a ~ 174° (slightly north of west), and the latitude of the Super-K detector 
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is 9 L ~ 36.3°, 

8m 



m l = — ^^°- 4eV2 - ( 2 °) 

For MINOS, a ~ 124° and L ~ 48°, so that 

ml ~ 0.008 eV 2 , (21) 

which is nearly two orders of magnitude smaller than the value required to describe the K2K data. 
The reason K2K gives a much smaller value for cos 2 (and hence requires a much larger value 
for ttiq) is that the neutrino path is nearly perpendicular to the Earth's rotation axis. Since the 
same m\ applies to both, £ = is excluded by a combination of the K2K and MINOS neutrino 
experiments. 

We note that for upward or downward atmospheric neutrino events, cos 2 = sin 2 6l, so m\ = 
5m 2 / sin 2 9l — 0.007 eV 2 , which is very close to the value extracted from the MINOS data. 

3.3 £^0 

If £ 0, then the preferred direction is not aligned with the Earth's rotation axis. For upward 
or downward atmospheric events there will be variation in cos 2 O (and hence in drn^j) over the 
sidereal period (see Eqs. (fT6j) and ([TT]) ). At the time of the sidereal day when <f> = or ir, there 
is always an extremum in cos 2 0. If |tan£| > |tan^|, then there are two more extrema at 
cos (j) = — tan#/,/ tan£. Thus there are two cases: 

• For |tan£| < |tan#£,|, the only extrema of cos 2 © occur at = and tt. Specifically, if 
< £ < 0l, then there is a minimum at (j) = and a maximum at <j> = tv, and 

sin 2 (# L -0 < c °s 2 © < sin 2 (# L + f ) . (22) 

If 7r — 6l < £ < ir, then the positions of the maximum and minimum reverse, and 

sin 2 (£ + 6 L ) < cos 2 < sin 2 (£ - L ) . (23) 

• For | tan £ | > |tan^| (i.e., 9l < £ < vr — 6£), cos 2 = when cos^ = — tan^/tan£ 
(which occurs twice a day) and there are maxima at </> = and tt with cos 2 = sin 2 (£ ± 6l). 
Therefore 

< cos 2 © < max[sin 2 (£-<9 L ),sin 2 (£ + f9 L )]. (24) 

The solid curves in Fig. [2] show the maximum and minimum values of cos 2 versus £ for upward 
and downward atmospheric neutrinos. For 0l < £ < vr — 9l, there are always two times during 
the sidereal day when cos 2 = 0, and hence there are no oscillations for up/down events (since 
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Figure 2: Maximum and minimum daily values for cos 2 G versus £ for MINOS (solid curves), K2K 
(dashed) and Super-K up/down atmospheric (dotted) data. In regions with two different local 
maxima, the larger one is shown. 

8 m eff = m o cos2 ©)■ This effect might be evident in the Super-K data if it were binned according 
to sidereal time. For values of £ less than 9l (or more than it — 6£), cos 2 G is always finite, with 
the degree of modulation decreasing as £ — > (or n). 

There is a similar situation for horizontal events, except that the critical angle that determines 
the number of extrema (and the values for the extrema) is 7 = sin -1 (sin a cos 9£) instead of 9^. 
For K2K, 7 ~ 5°, and the minimum cos 2 G is zero everywhere in the range 5° < £ < 175°. For 
< £ < 5° and 175° < £ < 180°, the minimum cos 2 Q is never larger than sin 2 5° ~ 0.008, so that 
there is always a time of day for K2K at which Sm^jj is suppressed and there are effectively no 
oscillations. The maximum and minimum cos 2 Q for K2K are also shown in Fig. [2j For MINOS, 
9l = 47.8° and the neutrino direction is approximately a = 124°; then 7 = 34° and the MINOS 
daily ranges for cos 2 are almost identical to those for Super-K up/down atmospheric events (see 
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Fig. ED. 

For £ < 7r/2, the maximum value for Sm^f in K2K is mgsin 2 ^ + 7^2^)) and for £ < jminos, 
the minimum value for Sm 2 ^ in MINOS is to 2 , sin 2 (£— Aminos)- For £ < 12° there is no value of m 2 , 
that gives Si^lff within both allowed experimental ranges (1.9 x l(r 3 eV 2 < 5m 2 < 3.5 x 10" 3 eV 2 
for K2K and 2.3 x 10~ 3 eV 2 < 5m 2 < 3.4 x 1CT 3 eV 2 for atmospheric neutrinos, at 90% C. L.)S 
Therefore, in an argument similar to the £ = case, the predicted MINOS and K2K firn 2 ^ disagree 
for £ < 12°, in contradiction with data, so that these values are excluded. For 12° < £ < 90°, there 
are always two times during the sidereal day when cos 2 O = for K2K, and there are no oscillations. 
For 12° < £ < 36°, atmospheric up/down events should show a significant modulation of 5m 2 ^, 
and for 36° < £ < 90° there are always two times during the sidereal day when cos 2 O = for 
atmospheric up/down events. Similar comments can be made for the range ir/2 < £ < ir. 

The results for K2K, MINOS and up/down atmospheric neutrinos may be summarized as fol- 
lows: 

• The range < £ < 12° (and by similar arguments, 168° < £ < 180°) is excluded by a 
comparison of the measured 5m 2 values in MINOS and K2K data. 

• For 12° < £ < 168°, there are always two times during the sidereal day when K2K should have 
no oscillations, i.e., no suppression of events relative to expectation. Up/down atmospheric 
neutrinos always have a significant modulation of 5m 2 ^, and for 36° < £ < 144° there are 
always two times during the sidereal day when up/down atmospheric neutrinos should also 
have no suppression. 

For horizontal atmospheric neutrino events {(3 = tt/2), cos© is given by Eq. (|18p : the daily 
fluctuations then depend on the compass direction of the event, a. Super-K has measured the 
compass dependence [10] and found agreement with an east-west asymmetry due to the Earth's 
magnetic field, plus oscillations. Any additional compass dependence must not be too large to 
remain consistent with the data. Table Q] shows cos 2 O for some typical values of a and (j). The 
direction dependence would not enhance or suppress the east-west difference, but could enhance 
or suppress oscillations along the east/west direction compared to north/south. Furthermore, 
enhancements could change to suppression (and vice versa) during the sidereal period. A detailed 
analysis would be needed to determine the compass-direction dependence for horizontal atmospheric 
neutrino events. 

1 This is approximately equal to the region where the cos 2 6 values do not overlap in Fig. [2] the difference is due 
to the slightly different ranges for Sm 2 in the two experiments. 
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Table 1: Values of cos 2 G for various values of a and (p. 



a 


= = 7r/2 = vr = 3vr/2 


0, 7T (E,W) 
vr/2, f (N,S) 


sin 2 £ sin 2 £ 
cos 2 (£ - 0l) cos 2 £ cos 2 0l cos 2 (£ + 6*l) cos 2 £ cos 2 6>l 



4 Solar neutrinos 

4.1 Directional dependence 

In a coordinate system (X', Y', Z') where the Z' axis is perpendicular to the Earth's orbital plane 
(the ecliptic plane), the direction of neutrino propagation may be written as (see Fig. [3]) 

p = cos tp X' + sin tp i , (25) 

where ip gives the position of the Earth in its orbit (ip = at the vernal equinox, tp = ir/2 at 
the summer solstice, etc.). The equatorial coordinates are related to the ecliptic coordinates via 
rotation by an angle r] ~ 23° about the X' axis, where r/ is the tilt of the Earth's rotation axis from 
the perpendicular to the ecliptic (see Fig. [3j). Then in the celestial equatorial frame the direction 
of propagation for solar neutrinos is 

p = cos tp X + sin ip cos r] Y — sin ip sin r/ Z , (26) 

and therefore 

cos = p ■ h = cos rp cos x sin £ + sin^(sinx sin£ cos r\ — cos£ sin 77) . (27) 

Note that cos 6 for solar neutrinos is independent of detector latitude (6l) and time of day (<p). 

4.2 Oscillation probability 

In matter there is an additional term in the hamiltonian due to coherent forward scattering of z/ e 's 
with electrons in matter, so that —2cE in the upper left element of h e ff is replaced by —2cE + 
V2GpN e and the mixing angle in Eq. ([6]) is then given by 

gzg£jjM 1 , (28) 

yj (cE - GpN e / y/2) 2 + a 2 cos 2 9_ 

where N e is the electron number density. For adiabatic propagation in the sun the solar neutrino 
oscillation probability is 

P(u e -► u e ) = cos 2 6 cos 2 Q + sin 2 6 sin 2 O , (29) 
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Figure 3: Position of the Earth in the ecliptic plane. The arrows represent the Earth's rotation axis 
and i/j = corresponds to the vernal equinox. The orientation of the celestial equatorial coordinate 
system relative to the ecliptic coodinate system is also shown. 

where 6q is the mixing angle at the creation point in the sun (with electron number density N® ~ 
90iV,4/cm 3 ) and 9 is the mixing angle in vacuum. P — > i at low energies! 2 ] There is a minimum in 



P at 



with minimum value 



E,, 



Pmin {ye 



2V2 c ' 
4a 2 cos 2 



< 



(30) 



(31) 



eJ 8a 2 cos 2 G + (G F N°) 2 2 ' 
where cos is given by Eq. f)2T[) . At E = 2E m i n there is a resonance and the probablity is \ , and 
for E > 2E m i n the probability increases monotonically, with limiting value unity as E — > 00. The 
angle depends on the time of year; averaging over t/j gives 



{Pmin {ye 



^{G F N»Y + 8a 2 D 2 



where 



D 2 = cos 2 x sin 2 £ + (sin x sin £ cos 77 — cos £ sin 77) 2 . 



(32) 



(33) 



If the probability minimum lies in the middle of the B solar neutrino region, then {P m in) m Eq. ([32 
will give the approximate survival probability of the 8 B neutrinos. 



2 The actual solar neutrino oscillation probability at low energies is closer to 0.7 |1X] . However, if an additional 
ajT e T term is included in h e f / , then the low energy probability can be fit to the higher value, at the expense of adding 
a fifth parameter to the model. 
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The formulas used above for the solar neutrino probability assumed adiabatic propagation. It 
can be shown that the propagation is adiabatic except close to the two times during the year where 
cos 6 = 0: 

. .x/cosxsinA . _ 1 /cos%sinA 
tp = — sin I — I and ir — sin I — J ; (34) 

this was also pointed out in Ref. [5] for the special case £ = 0. To include the effects of nonadiabatic 
propagation, Eq. (f29l) must be modified to 



P(v e i/ e ) = - [l + (1 - 2P x )(cos 2 6 cos 2 9 + sin 2 sin 2 6 Q )] , (35) 

where P x is the level-crossing transition probability, 

P x = ( (36) 

and 7 r is the adiabaticity of the transition at the level-crossing resonance. For our Hamiltonian 

_ 2V2a 2 cos 2 e 
7r G F \dN e /dL\ r ' 1 j 

where \dN e /dL\ r is the rate of change of iV e at the resonance. At E m i n the probability becomes 

mml e ej " 8a 2 cos 2 6 + (G F N»Y + (GpN®) 2 + 8a 2 cos 2 6 x ' {6 *> 

where the first term on the right-hand side is the adiabatic contribution and the second term the 
nonadiabatic correction. Propagation is nonadiabatic when 7 r is small, which occurs when a 2 cos 2 G 
is small. For the parameter ranges of interest we find that 8a 2 cos 2 O <C (GfN®) 2 in the regions 
where P x is nonnegligible, so that the probability reduces to 



From Eqs. (|36|) . (|37j) and (|39j) we see that the survival probability goes to unity when cos© = 0. 
The probability is shown in Fig. Inversus V using both the adiabatic and nonadiabatic formulas; 
they differ substantially only near the values of if) given by Eq. (|34p . 

4.3 Constraints from solar data 

In order to fit the solar neutrino data, (Pmin) must match the measured probability for the 8 B 
neutrinos, i.e., (P m in) — 0.34 (we use the ratio of CC to NC rates in SNO [12] to avoid complications 
due to theoretical uncertainties in the solar neutrino spectrum). Since there is no apparent energy 
dependence in the 8 B oscillation probability, the minimum must occur near the middle of the 8 B 
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Figure 4: Representative solar neutrino survival probability at E m in using the adiabatic (Eq. (|3lj) . 
solid curve) and nonadiabatic (Eq. dashed) formulas, shown versus the time of year (measured 
by tp). The parameters for this example are a = 7 x 1CP 12 eV, £ = 45° and x = 0, which give 
(Pmin) — 0.34. The two formulas differ only close to the values of ip where cos = 0. 

spectrum (E m in — 10 MeV), so that probabilities at either end of the spectrum are not much larger 
than in the middle. This results in the two constraints (from Eqs. (|30p and (I32jl ) 

G F N® 



2\/2 Emin 



1.7 x 10 



-19 



iD ~ V2G F N 



V / (^"i*«)(l (Pmin)) _ g q x ]^Q-12 g y 



(40) 
(41) 



(1 - 2(P min )) 

where Eq. (|41|) uses the adiabatic expression for (P m in)- Since Eq. ([ID]) depends only on the initial 
density for 8 B neutrinos and the central energy of the SNO spectrum, we will use this result for c 
throughout the rest of this paper. 

We note that although the value of a required to fit (Pmin) depends on the value of D (which in 
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turn depends on the preferred-direction parameters £ and x), the product aD is fixed by Eq. (|41j) . 
and the oscillation probability versus time will always be identical to that shown in Fig. [4j except 
for a possible shift in phase and the corrections for the two nonadiabatic spikes. This can be 
understood by rewriting the adiabatic probability in Eq. (|3ip as 

4a 2 D 2 sm 2 (ip + 5) 



8a 2 D 2 sin 2 (V + 5) + (G F iV e ) 2 

where 



(42) 



tan 5 ^ Sil ^ C ° SX . (43) 

sin x sin £ cos r\ — cos £ sin r\ 

Thus the time variation of P m in has the same shape and maximum and minimum values when aD 
is held fixed. 

The measured solar neutrino survival probability does not exhibit much variation throughout 
the year. The SNO collaboration has tested their solar neutrino data for periodicities [13] and 
found a variation during the year that is consistent with the 1/r 2 dependence of the flux as the 
Earth's distance from the sun varies. The uncertainties in the rate are of order 3-5%, so there is 
little room for any additional annual variation. The SNO periodicity data sample includes all of 
their solar neutrino data in both the D2O phase and salt phase, and combines events from charge- 
current (CC), neutral-current (NC), electron scattering (ES) and backgrounds (B). They measured 
the relative event rate versus time of year, normalized to the mean rate, i.e., 

r> = Kc + 3cL + n es\p + *■(! - n + gg (u) 

N» NC + N° c (P)+N° ES [(P) + r(l-(P))]+N°' 1 ' 

where P is the oscillation probability, N® is the number of events expected without oscillations, r 
is the ratio of the NC to CC cross sections and angle brackets indicate mean values. For {P m in) = 
0.34, the bicycle model with directional dependence predicts R should vary between 0.42 and 
1.19 throughout the year. Since the SNO measurement of R varies by at most 5% at any time 
during the year, the pure direction-dependent case clearly cannot fit the SNO periodicity test while 
simultaneously reproducing the correct average survival probabilty. 

To verify this quantitatively we have searched the a, £ and x parameter space via Monte Carlo, 
using the twelve bins of the SNO periodicity data and the SNO average probability (0.34 ± 0.03, 
from the CC to NC ratio). We have used the appropriate weighting of run times and D20/salt 
phase for each bin, and used Eq. (|39|) for the oscillation probablity, which includes the nonadiabatic 
part. The nonadiabatic spikes appreciably affect bin-averaged probabilities only in the bins where 
they occur, and then by order 0.05 or less. In Fig. [5] we show the SNO periodicity data plus the 
best fit when varying over a, £ and x when (P) is constrained to lie within la of the central value. 
Allowing (P) to lie outside the la range can improve the x 2 > but in all cases the x 2 per degree 
of freedom (DOF) is such that the probablity that the model describes the data is 2 x 10 -8 or 
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Figure 5: Best-fit prediction for R in Eq. (|44h for the 8 B neutrinos (solid curve) and the SNO 
measured value for R (data points) shown versus time of year (measured by ip). Both the SNO 
data and the model predictions are averaged over each bin, and the SNO data has been corrected 
for the 1/r 2 variation due to the changing Earth-Sun distance. The prediction for (P) has also 
been constrained to lie witin la of the SNO central value. The model parameters for the best fit 
are a = 1.96 x 10~ 12 eV 2 , f = 43° and x = 298°, with x 2 /DOF = 361/10. 



less. The best fit has very little annual variation, but (P) ~ 0.21. Therefore we conclude that 
the generalized direction-dependent bicycle model is strongly ruled out solely by the solar neutrino 
data. 
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5 Combined constraints 



5.1 Adding a direction independent term 

Since the pure direction-dependent case is ruled out, we now generalize the model to include both 
direction-independent as well as direction-dependent terms in the off-diagonal elements of h e ff. 
This increases the number of parameters in the model to five. If we define (a^,)^ = (oi)c T = 
(a cos p, a sin p n/y/2), where n is again the preferred direction, then cosG should be replaced by 
cos p + sin p cos in our previous formulas. The parameter p determines the amount of direction 
dependence: /9 = 7r/2or37r/2 corresponds to the pure direction-dependent case we discussed before, 
while p = or n corresponds to no direction dependence. 

For a given preferred direction (fixed £ and x), the parameters c and a are determined from 
the solar neutrino data using Eqs. ([30]) and ([3"T]) . after the substitution cos© — ► cosp + sin p cos 
is made. Then using Eqs. ([5]) and (fTTj) . 5m 2 , for long-baseline and atmospheric neutrinos may be 
written as 

a 2 

5m 2 ^ = — (cos p + sinpcos ©) 2 . (45) 



It is convenient to rewrite cos as 

cos G = D sin(V> + 6) , (46) 



where 5 is defined in Eq. (|4"3"]l. Integrating P m in in the modified Eq. (pTTj) over ip, leads to 

X 1 mini — 2 



1 - ^0t\] s2 + 8a2 ( D2 sin2 P ~ cos2 P) + {G F N2) 2 



for adiabatic neutrinos, where 



(47) 



S 2 = y 64a 4 (L> 2 sin 2 p - cos 2 p) 2 + 16a 2 (G F iV e ) 2 (D 2 sin 2 p + cos 2 p) + (G F iV e ) 4 . (48) 
5.2 No direction dependence 

For the pure direction-independent case (p = or tt), 5m 2 ^ = m 2 , = a 2 /c for atmospheric and 
long-baseline neutrinos and P m in for 8 B solar neutrinos is given simply by Eq. (|3ip ; for P m %n = 0.34, 



Pmin G F N® | 2 



2(1 - 2P„ 



2.5 x 10~ i2 eV , (49) 



and the prediction from the solar neutrino data is Sm 2 ^ = 3.6 x 10 -5 eV 2 for atmospheric and long- 
baseline neutrinos, which is clearly in contradiction with the data. Therefore, the pure direction- 
independent case is ruled out by the combined data. 
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5.3 Mixed case 

For a mixture of direction-dependent and direction-independent terms in h e ff, a fit must be done 
to the solar data to determine an allowed region in parameter space, and then the predictions 
for dm^ff in long-baseline experiments can be compared to data. To fit the solar data we take 
the 12 bins from the SNO periodicity data sample for the relative rate R and add the additional 
constraint that the average oscillation probability must be P = 0.34 ± 0.03, as described in the 
previous section. As before, we fix the value of c to that given in Eq. (|40p . and vary over the 
parameters £, %i P an d a with a Monte Carlo. The 99% C. L. allowed regions are determined by 
restricting the x 1 /DOF to be less than 2.4 for nine DOF (there are thirteen data points and four 
parameters). The best fit to the SNO data has a = 3.0 x 10~ 12 eV 2 , £ = 21°, % = 94° and p = 114°, 
with x 2 /DOF = 4.84/9. 

Predictions for 5to 2 jj can then be made for K2K and MINOS. Since 5m 2 ^ depends on cos O, 
it will vary during the sidereal day for £ ^ 0, with ranges depending on £ as shown in Fig. [2j The 
strictest constraints come from K2K; maximum possible values of <5 m e// ^ n K2K are shown versus 
£ in Fig. El 

In most all cases the maximum possible dm^f can never be in the experimentally measured 
range 1.9 x 10~ 3 < Sm 2 < 3.5 x 10 -3 eV 2 . Only a small region near £ ~ 7] = 23° or it — rj ~ 157° 
can give a large enough value of <5m 2 ^. This allowed region is also characterized by p ~ tt/2 or 
37r/2 and x — 7r /2 or 3ir/2, such that |-Dsin/?| <C | cosp\, and values of a > 3 x 10~ n eV. 

As evident from Eq. (|46|) . this results in cos ~ for solar neutrinos (i.e., the preferred direction 
is nearly perpendicular to the ecliptic plane), so that the directional dependence for solar neutrinos 
is minimal, even though the direction-dependent coefficient sin p is much larger than the direction- 
independent coefficient cos p. For atmospheric and long-baseline neutrinos this fortuitous situation 
does not occur and the direction-dependent piece is sizable, with daily variations of cos 2 O given by 
Fig. [5J Therefore the case with a mixture of direction-dependent and direction-independent terms 
is severely constrained, and there is a strong variation of Sm 2 ^ for atmospheric and long-baseline 
neutrinos during the sidereal day for the allowed solutions. 

5.4 KamLAND 

For reactor neutrinos with both direction-dependent and independent terms, from Eqs. dl}-® we 
have 

P(v e v e ) = 1 - sin 2 29 sin 2 Q A 31 Lj , (50) 

where 

A 3 i = 2 v / (c£) 2 + a 2 (cosp + sinpcos6) 2 , (51) 
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Figure 6: Maximum value of 5m 2 in K2K allowed at 99% C. L. by model parameters consistent 
with the solar R and (P) data, shown versus the preferred direction £. 

sin 2 29 = l-T—yi 37 ; ^ , (52) 

[chjy + a z (cos p + sin p cos Hj z 

and cos is given by Eq. Q18p . For the values of the parameters that fit solar data and give a large 
enough 5m 2 , for long-baseline neutrinos, (cE) 2 <C a 2 and sin 2 29 ~ 1 at reactor neutrino energies, 
except possibly for the brief time of day when cos p + sin/? cos 0^0. 

As discussed in Sec. 3, given 7 = sin - 1 (sin a cos 9l), where a is the compass direction of the 
incoming neutrino and 9l the latitude of the detector, the maximum and minimum values for 
cos 2 during the sidereal day are given by Eqs. (|22|) - (|24|) with 9l replaced by 7. Then for the 
parameter ranges found above, it is not hard to show that for all of the reactors contributing to 
the KamLAND signal, the oscillation argument varies over many cycles during the day, so that 
the oscillation probability is close to 0.5, regardless of neutrino energy. Thus the bicycle model 
gives a suppression in KamLAND that is nearly independent of energy, contrary to the KamLAND 
data [H], which excludes an energy-independent suppression at 99.6% C. L. 
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We have verified this result numerically using typical solar /atmosperic/long-baseline solutions, 
averaging over the sidereal day, and summing over individual reactor contributions - the suppression 
varies by at most 0.02 over the range 2.5 < E < 6 MeV which supplies the bulk of the KamLAND 
data. The average survival probability of the bicycle model solutions is at most about 0.55, well 
below the measured KamLAND value of P = 0.658 ± 0.044 ± 0.047. Therefore the bicycle model 
with a mixture of direction-dependent and direction-independent terms is also excluded. 

6 Conclusions 

We have shown that the generalized five-parameter bicycle model with Lorentz-invariance violation 
and no neutrino masses can be ruled out by a combination of solar, long-baseline and reactor 
neutrino data. The pure direction-dependent case is ruled out because it gives a large annual 
variation in the oscillation probability for 8 B solar neutrinos, at odds with SNO periodicity data. 
The pure direction-independent case is ruled out because the values of the parameters required to 
fit the SNO data predict a value of 5m 2 in long-baseline experiments that is too small by nearly two 
orders of magnitude. Having a mixture of direction-dependent and direction-independent terms in 
the off-diagonal elements of h e ff is excluded when KamLAND is added to a combination of solar 
and long-baseline data. 

Although the five-parameter bicycle model cannot fit all of the data, the full h e ff with Lorentz- 
noninvariant oscillations of massless neutrinos has 160 parameters [5J, and a comprehensive com- 
parison with data is impractical. However, it is clear that any direction dependence will encounter 
severe constraints, including variations during the sidereal day which were not pursued in this paper. 
Restricting h e ff to only direction- independent terms reduces the number of Lorentz-noninvariant 
parameters to 16 [5]. Even then, as our analysis of the direction-independent bicycle model sug- 
gests, finding a set of parameters that would simultaneously fit solar, atmospheric, long-baseline 
and reactor data will be difficult at best. 
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